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Abstract The solar atmosphere is magneticaUy structured and highly dynamic. 
Owing to the dynamic nature of the regions in which the magnetic structures 
exist, waves can be excited in them. Numerical investigations of wave propa- 
gation in small-scale magnetic flux concentrations in the magnetic network on 
the Sun have shown that the nature of the excited modes depends on the value 
of plasma /3 (the ratio of gas to magnetic pressure) where the driving motion 
occurs. Considering that these waves should give rise to observable characteristic 
signatures, we have attempted a study of synthesized emergent spectra from 
numerical simulations of magneto-acoustic wave propagation. We find that the 
signatures of wave propagation in a magnetic element can be detected when the 
spatial resolution is sufficiently high to clearly resolve it, enabling observations in 
different regions within the flux concentration. The possibility to probe various 
lines of sight around the flux concentration bears the potential to reveal different 
modes of the magnetohydrodynamic waves and mode conversion. We highlight 
the feasibility of using the Stokes-F asymmetries as a diagnostic tool to study 
the wave propagation within magnetic flux concentrations. These quantities can 
possibly be compared with existing and new observations in order to place 
constraints on different wave excitation mechanisms. 

Keywords: Magnetic fields. Photosphere; Magnetic fields. Models; Magneto- 
hydrodynamics; Spectral Line, Intensity and Diagnostics; Polarization, Optical; 
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1. Introduction 

Spectral lines inform us on the properties of the atmosphere in which they form. 
In addition, any presence of magnetic fields in the atmosphere modifies the 
polarization state of the light emerging from the surface. A wealth of informa- 
tion about the structure and dynamics of the magnetized regions of the Sun is 
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hidden in the polarization state of spectral lines. In this paper, we will study the 
polarization signatures of magneto- acoustic wave propagation in a photospheric 
magnetic flux concentration using magnetohydrodynamic numerical simulations 
thereof. The initial and simulated models are thought to represent a magnetic 
flux concentration in the magnetic network of the Sun. 

The two-dimensional simulations by Steiner et al. (1998) and Grossmann- 
Doerth, Schiissler, and Steiner (1998) showed that the Stokes proflles vary strong- 
ly in response to a dynamic magnetic atmosphere. Using a similar forward 
modelling, the effect of waves on Stokes-F profiles of Ca ii infrared lines were 
studied by Pietarila et al. (2006). They saw a clear time-dependent behaviour 
of the Stokes-y proflles as a result of wave propagation and shock formation 
occurring in the numerically simulated atmosphere. Even though they were able 
to reproduce the atmospheric dynamics in the form of observational signatures in 
the Stokes profile, their work was limited to the weak field case. Several authors 
(Rosenthal et al., 2002; Bogdan et al., 2003; Cranmer and van Ballcgooijcn, 
2005; Hasan et al., 2005; Hasan and van Ballegooijen, 2008; Khomenko, Collados, 
and Felipe, 2008; Vigeesh, Hasan, and Steiner, 2009; Fedun, Erdelyi, and Shelyag, 
2009; Kato et al., 2011) have investigated wave phenomena in magnetically 
structured atmospheres. Shelyag et al. (2010) constructed a photospheric bright 
point model and studied the observational signatures of wave propagation in 
them including the response in Stokes V . The effect of a direct excitation of a 
magnetic flux concentration by granular buffeting and the corresponding spectral 
signatures of wave propagation and mode coupling and mode conversion in these 
structures have not been studied so far. Fujimura and Tsuneta (2009) report 
on the observation of magnetohydrodynamic wave propagating along magnetic 
flux tubes in the solar atmosphere. Their work is based on observations with 
the SOT/SP instrument onboard the Hinode spacecraft and demonstrates the 
feasibility of such measurements. 

Inspired by these observations and by numerical simulations, and driven by 
the desire to find possible polarimetric signatures of wave propagation in mag- 
netic network elements, we have attempted to study the feasibility of using the 
Stokes- spectra obtained from our simulation as a diagnostic for magnetohy- 
drodynamic wave propagation. 

The paper is structured as follows: Section 2 will give the construction of the 
initial equilibrium model. In Section 3, we will explain our numerical method and 
boundary conditions. In Section 4, we will present three different experiments 
of wave propagation in magnetic elements and in Section 5 we will discuss the 
properties of the Stokes proflles emerging from the simulation box. The summary 
and conclusions will be given in Section 6. 



2. Initial Equilibrium Model 

We construct a two-dimensional initial atmosphere in Cartesian coordinates 
containing a magnetic flux sheet. For the construction we use the numerical 
methods described in Steiner, Pneuman, and Stcnflo (1986). The magnetic field 
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configuration and the pressure distribution in the physical domain is specified 
as in Vigccsh, Hasan, and Steiner (2009). 

The magnetic field can be written in terms of the flux function Tp(x, z) as 



Bx — — 



chp_ 
dz ' 



B, = 



dx 



(1) 



Contours of constant flux value, i), correspond then to magnetic field lines. We 
identify V = with the symmetry axis in the centre of the fiux sheet and 
ij) — ±i/)iuax with the side boundaries, which also defines the total magnetic flux. 
The gas pressure is prescribed as a function of height and fleld line, p{ip,z), in 
the following way. 
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where the constants a, b, c, and d are chosen such that the pressure and its 
first derivative with respect to ip are continuous functions of ip and where we 
choose n = 8. Bq and po are the magnetic field strength and the gas pressure, 
respectively, on the axis of the flux sheet at the reference height z = 0. P2, tpi, 
and tjj2 are chosen conveniently so as to obtain the desired cross-sectional profile 
for the magnetic field, viz., Bz{x) at z = 0, as shown further below in Figure 3. 
The gas pressure along the axis is defined as. 



p(0, z) = Po exp \ - 



RT{z) 



(3) 



Here /i is the mean molecular weight, taken to be /Lt = 1.297, R the universal 
gas constant, and g the gravitational acceleration at the solar surface. T{z), the 
temperature as a function of height, is modelled by an analytical function of the 
form 



T{z) — Tq + atanh(7z + c). 



(4) 



It is constant on levels of constant z, assuming that an efficient radiative ex- 
change between fiux sheet and the ambient medium levels the temperature to 
a horizontally isothermal state. With an appropriate choice of Tq, a, 7, and 
c, we construct a photospheric temperature run as shown in Figure 1. For 
comparison. Figure 1 also shows the temperature as a function of height of 
the solar atmospheric model of Holweger and Mueller (1974), which was derived 
from observations under the assumption of local thermodynamic equilibrium 
(LTE). Initially, the temperature drops rapidly from 10 500 K at the bottom 
boundary to 6 300 K at z = km, then asymptotically decreases to 4000 K. 
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Figure 1. Temperature as a function of height of the model atmosphere, according to Equa- 
tion (4) using To = 9000 K, a = -5000 K, 7 = 3 X 10"* cm"!, and c = 0.6. z = corresponds 
approximately to continuum optical depth unity, Tc = 1, for 1/) = Vmax. The dashed curve 
shows the photospheric reference model of Holweger and Mueller (1974). 




This temperature profile approximately reflects the temperatme profile of the 
Holweger and Miiller model in the photospheric part and enables us to compute 
spectral lines in absorption. We do not include a chromospheric temperature rise 
because we intend to compute the spectral lines in LTE only. With the choice of 
an analytical expression for the temperature and an ideal gas equation of state 
with constant molecular weight, the initial atmosphere quickly relaxes to the 
static initial solution of the discrete numerical scheme, although at the expense 
of a more realistic atmosphere. 

Having defined the gas pressure and the temperature distribution through 
Equations (2)-(4), we obtain the density distribution. From the force balance 
perpendicular to the direction of the field lines, one obtains the electric current 
density, 



dp 



(5) 



The new magnetic field configuration can be calculated from the current density 
using the Grad-Shafranov equation, 

d^+9^= (6) 

A detailed derivation of these equations is given in Steiner (2007). 

Due to the symmetry of the problem, we solve Equation (6) in a computational 
domain that consists of only half of the fiux sheet of horizontal and vertical 
extensions of 640 km and 1 600 km, respectively. The bottom boundary is at a 
depth of z = —300 km. The domain is discretized on an equidistant rectangular 
mesh with a spacing of 5 km. The left side of the domain corresponds to the axis 
of the fiux sheet. The value of ijj is prescribed at the left and the at right side 
boundaries. At the top and bottom boundaries, we use the Neumann condition 
dip /dz = 0, assuming that the horizontal field component vanishes at these two 
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Figure 2. Magnetic field strength as a function of height on the axis (soHd curve) and in the 
ambient medium (dashed curve) of photospheric magnetic flux sheets with field strengths of 
1000 G and 1600 G at the axis at 2 = 0. 
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Figure 3. (top row) and (bottom row) components of the magnetic field as functions 
of horizontal distance at the following heights; z = km (solid curve), z = 500 km (dotted 
curve), and z = 1000 km (dashed curve). The plots refer to photospheric magnetic flux sheets 
with field strengths of (a) 1000 G (gauss) and (b) 1600 G at the axis at 2 = 0. 
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Table 1. Equilibrium model characteristics for the 1000 G and 1600 G flux sheets. 
The numbers in the first row of each physical quantity correspond to the top boundary 
{z = 1300 km) and the numbers in the second row correspond to the height z = km. 
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boundaries. Starting from a reasonable initial field configurations, e.g., as derived 
from the thin flux-tube approximation, one obtains by iteration of Equations (1)- 
(6), a final, sclf-consistcnt, magnctohydrostatic model. 

We calculate cases corresponding to different field strengths (at z = 0) ranging 
between 1000 G (gauss) and 1600 G, on the axis of the sheet. For the case with 
weaker field strength, the /3 = I layer lies well above the bottom boundary 
dropping to a minimum height of 370 km only. Here /3 is the ratio of the gas 
pressure to the magnetic pressure. Any field line originating from the bottom 
eventually crosses this layer, dividing the flux sheet into two regions; a lower 
region with /3 > 1 and the upper region with (3 < 1. We say that the flux sheet is 
rooted in a high (3 region. In the case of a stronger field, e.g., when Bq = 1600 G, 
the = 1 layer traces the boundary of the flux sheet. The entire flux sheet is in 
a region of /3 < 1, and only the ambient medium has a high /3 plasma value. 

The variation of the magnetic fleld strength with height on the axis and in 
the ambient medium is shown in Figiuc 2. In both cases, the magnetic fleld 
drops to a uniform value within the flux sheet. In the ambient medium, the 
fleld strength in the lower part is negligible, but with increasing height, it settles 
down to the same uniform value as that of the flux sheet. Figure 3 shows the 
horizontal variation of the horizontal and vertical components of the magnetic 
fleld at three different heights, z = km (solid curve), z = 500 km (dotted 
curve), and z = 1000 km (dashed curve). The horizontal component of the field 
at z = km is close to zero, hence the fleld is almost vertical at this level. The 
flux sheet at this height has a vertical component of the magnetic field that 
drops sharply to the ambient value in the horizontal direction, confining it to a 
narrow region with a width of about 320 km. The fiux sheet expands with height 
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to cover the entire horizontal extent with a homogeneous vertical field, starting 
at a height of below z — 1000 km. 

The equilibrium characteristics of the two models are summarized in Table 1. 
The values of the physical quantities are given at the sheet axis and in the 
ambient medium for the top boundary (z = 1300 km) and for the height z — km 
in the first and second rows of each entry, respectively. Note that the plasma /3 
at the sheet axis at the base is 2.1 for the 1000 G case and 0.2 for the 1600 G 
case. The sound speed (c^) and the Alfven speed {va) are defined as, 




B 

- (8) 



where 7 is the ratio of specific heats taken to be 5/3 and p, p, and B are 
the equilibrium values of gas pressure, density and magnetic field strength, 
respectively. 



3. Numerical Simulation: Methods and Boundary Condition 

Wave propagation is studied by an impulsive transverse excitation of the lower 
boundary in the equilibrium model (similar to Hasan et al. (2005) and Vigccsh, 
Hasan, and Steiner (2009)). The system of MHD equations, given in conservation-l 
law form for an inviscid adiabatic fluid, is solved according to the method 
described in Steiner, Knolker, and Schiissler (1994). 

The side boundaries are open due to a constant extrapolation of the variables 
from the physical domain to the boundary cells. The horizontal component of 
the momentum at the top and bottom boundaries and the vertical component 
at the top boundary are also set by a constant extrapolation. The density in the 
top and bottom boundary cells is determined using a linear-log extrapolation. 
For the temperature, constant extrapolation is used at the top boundary. The 
temperature in the bottom boundary cells is determined using Equation (4) . The 
horizontal component of the magnetic field at the top and bottom boundaries 
are set equal to the corresponding values at the preceding interior point so that 
dBx/dz — 0. The vertical component of the magnetic field is determined by 
the condition V • B = 0. With these boundary conditions, acoustic waves (slow 
modes) can propagate across the top boundary with little reflection. They also 
work for the fast, magnetically dominated mode as it typically gets refracted and 
converted to an acoustic mode before leaving the computational domain across 
the side boundaries as will be shown in Section 4. Thus, only a small fraction 
of the fast-mode wave fronts reaches the top boundary: the bulk of it leaves the 
computational domain through the side boundaries. 
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Similar to Vigccsh, Hasan, and Steiner (2009), the transverse velocity Vx at 
z = is specified as follows: 



where Vq denotes the amplitude of the horizontal motion and P is the wave 
period. This form simulates an impulsive transverse excitation of the flux sheet 
at the lower boundary. For simplicity, we assume that all points of the lower 
boundary have this motion: this does not generate any waves in the ambient 
medium, other than at the interface with the flux sheet. In order to achieve 
significant intensity signals, we use Vg = 5 kms~^ and P = 24 s. Such short 
duration motions are expected to be generated by the turbulent motion in the 
convectively unstable subsurface layers, where the flux sheet is rooted. Cranmer 
and van Ballcgooijen (2005) studied the kinematics of G-band bright points 
and suggested that there are two components involved: a "random walk phase" 
and a "jump phase" . Our work considers the case with higher velocities, which 
represents the "jump phase" component of Cranmer and van Ballegooijcn (2005). 
This motion generates magnetoacoustic waves in the flux sheet. 

4. Dynamics 

We consider a unidirectional horizontal displacement of the entire bottom region 
below z = km (wide excitation region), first for the case in which the field 
strength is 1000 G at z = (moderate field case) and second, for the case in 
which the field strength is 1600 G at z = (strong field case). For the strong 
field case, we also consider a narrow excitation region with excitation below 
z = —150 km only. The excitations correspond to the impulsive case discussed 
in Section 3 and given by Equation (9). 

4.1. Moderate Field Case 

The /3 = 1 contour in this case is well above z = and hence all the magnetic 
field lines that emerge from the base cross this layer at some height as is visible 
in Figure 4. 

The motion in the direction to the right hand side takes place in the region 
where /3 > 1 (high-/? excitation) and results in the excitation of waves in the form 
of a fast (predominantly acoustic) wave and a slow (predominantly magnetic) 
wave, which propagate at the sound speed and the Alfven speed, respectively. 
The fast wave manifests itself as a compression and rarefaction of the gas at 
the leading and trailing edges of the flux sheet, respectively. It can be clearly 
discerned in the snapshot of the temperature perturbation, ST (the temperature 
difference with respect to the initial value), shown in Figure 4, which is taken 
60 s after the start of the perturbation. The black curves denote the magnetic 
field lines and the white curve depicts the /3 = 1 contour. The perturbations 
are 180° out of phase on the opposite sides of the sheet axis, leading to a quasi 
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Figure 4. Temperature perturbations and velocity field of a flux sheet in which the field 
strength at the axis at z = is 1000 G. The snapshot is taken at time 60 s after the initiation of 
an impulsive horizontal motion of the entire region below 2 = km according to Equation (9). 
The amplitude of the motion is 5 kms~^ and the period P = 24 s. The black curves represent 
magnetic field lines and the white curve depicts the /3 = 1 contour. Vertical lines indicate 
different bundles of lines of sight considered for the Stokes analysis. 

anti-symmetric wave pattern. As these fast waves travel upwards they eventuahy 
cross the layer of /3 = 1, where they change from fast to slow, without changing 
their acoustic nature: this corresponds to a "mode transmission" in the sense 
of Cally (2007). The transmission coefficient depends (among others) on the 
"attack angle", i.e., the angle between the wave vector and the local direction of 
the magnetic field (Cally, 2007). On the (3 = 1 layer, away from the sheet axis, 
where the wave vector is not exactly parallel to the magnetic field, we do not 
have complete transmission of the fast wave to a slow wave. Rather, there is a 
partial conversion of the mode from fast acoustic to fast magnetic, so that the 
energy in the acoustic mode is reduced correspondingly. For more details, the 
reader is referred to Vigeesh, Hasan, and Steiner (2009). 

4.2. Strong Field Case 

We consider a uniform horizontal displacement of the entire bottom boundary 
region of a thickness of 150 km (which we term as narrow) and of a thickness of 
300 km (wide) from the bottom boundary. This is to mimic the buffeting of the 
flux tubes by granular eddies of two different depths. The excitation corresponds 
to the impulsive case as given by Equation (9). Figure 5 shows the temperature 
perturbation 6T at 40, 60, and 80 s for the two cases of narrow and wide regions 
of excitation. 

Here, the contour of /? = 1 approximately traces the boundary of the flux 
sheet (see Figure 5). The transverse motion of the lower boundary generates 
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Figure 5. Temperature perturbation for a narrow and a wide excitation depths for a flux 
sheet in which the field strength at the axis at 2 = is 1600 G. Time instances refer to 40, 
60, and 80 s (from bottom to top) after the initiation of an impulsive horizontal motion in the 
regions from (a) z = —150 km and (b) 2: = km to the bottom boundary. The duration of the 
motion is 12 s and has an amplitude of 750 ms~^. The thin, black, vertically running curves 
represent field lines and the white curve corresponds to the contour of /3 = 1. The arrows 
indicate velocities of 50 ms~^ and more. Vertical lines indicate different bundles of lines of 
sight considered for the Stokes analysis. 



slow (predominantly acoustic) and fast (predominantly magnetic) waves. Since 
the contour of /3 = 1 runs along the boundary of the flux sheet, the waves 
that travel within the flux sheet along the magnetic field lines upwards do not 
encounter this layer and hence do not undergo mode conversion. On the other 
hand, the fast wave, which can travel across the field, encounters the /3 = f 
contour at the boundary of the flux sheet. As the fast wave crosses this boundary. 
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it enters a region of negligible field and hence gets converted into a fast (acoustic) 
wave as can be seen in the snapshot of temperature perturbations at an elapsed 
time of 40 s. At this time, the fast wave in the low-/3 region, which is essentially a 
magnetic wave, undergoes mode conversion and becomes an acoustic wave, which 
creates fluctuations in temperature visible as wing-like features at the periphery 
of the flux sheet (approximately along the (3 = 1 contour). The refraction of 
the fast wave, due to the gradients in the Alfven speed, and the eventual mode 
conversion are prominent in the case of the wide excitation regime. It can be 
hardly seen in the case where the flux sheet is shaken over a narrow region. 
This is due to the fact that a wide excitation range in a flux tube embedded in a 
low-/? region creates more magnetic pressure fluctuations relative to gas pressure 
fluctuations and hence imparts more energy to the fast (magnetic) wave. Due to 
the gradients in the Alfven speed, this mode gets refracted and returns back to 
a region with high-/?, when this energy is transferred to the fast (acoustic) wave, 
producing larger temperature fluctuations. This has implications for a realistic 
atmosphere, where granular eddies of different sizes are likely to impact deep 
rooted flux tubes. An impact over a wide range on a flux tube with strong field 
will transfer more energy to the fast (magnetic) mode. The ambient atmosphere 
regains part of this energy in the form of a fast (acoustic) wave due to the 
refraction of the fast mode and eventual mode conversion. When the excitation 
range is narrow, there is relatively weak magnetic pressure fluctuation compared 
to gas pressure fluctuation and hence, most of the energy goes into the slow 
(acoustic) mode, which is channelled up along the flux tube and eventually 
dissipates by shock formation. 

5. Stokes Diagnostics 

Even the largest modern solar telescopes are still not capable of resolving small- 
scale magnetic structures completely. Observations of the Stokes parameter / 
usually cannot reveal properties exclusive to the magnetic feature, since they 
suffer from mixing contributions from the magnetic field with contributions from 
the surrounding field-free or weak- field plasma. Differently from that, Stokes V 
gives us the properties intrinsic to the magnetic structure, because the circularly 
polarized light is formed only where the magnetic field is present. Therefore, 
presently the most sensitive method to study the magnetic atmosphere is by 
analyzing the Stokes-^ spectra emerging from them (see Sigwartli (2000) for a 
review) . 

We have computed the emergent Stokes-V^ profiles from the top of our simula- 
tion box for the weak field and the strong field cases, using the Stokes radiative 
transfer code DIAMAG (Grossmann-Doerth, 1994). This code calculates the 
normalized Stokes parameters by solving the Unno-Rachkovsky equations of 
radiative transfer. At the same time, it computes the line depression contribution 
function for each wavelength point. The program requires the temperature, gas 
pressure, magnetic field vector, velocity, and micro-turbulence to be specified 
on every grid point along the line of sight. The calculations were done for a 
set of four Fe i lines, viz., XX 5250.2, 5247.05, 6301.5, and 6302.5 A. The atomic 
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Table 2. Atomic parameters of the selected lines. 
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Values taken from Nave et al. (1994). 



parameters of the selected lines are listed in Table 2. We set the microturbulence 
velocity to zero for all the line-transfer calculations. 

Table 2 lists the wavelengths of the lines, the excitation potentials of the lower 
level, oscillator strengths (log ((?/)), and the effective Lande factors (^eff)- The 
two pairs of lines were selected because each forms under similar conditions in 
the atmosphere, since the lines of each pair have similar excitation potentials 
and oscillator strengths, which means similar opacities. But the difference in 
Lande factor for these lines make them useful for measuring the magnetic field 
strength, in particular, deviations from the weak field regime. These lines are 
commonly used to study solar magnetic fields. Socas-Navarro et al. (2008) have 
confirmed the reliability of using these four lines for the diagnostics of the quiet 
Sun magnetic field. 

Here we study the spectral signature of wave propagation in magnetically 
structured atmospheres with dynamically varying magnetic field. The effects of 
wave propagation in the four Fe i lines, listed in Table 2 are assessed, using 
the numerical simulations described in Section 4. The Stokes spectra for vertical 
lines of sight, separated by a horizontal distance of 10 km, were computed for 
each time step. These correspond to real observations at disk centre. Here, we 
present the analysis of the Stokes-V^ spectra for the two cases of moderate and 
strong magnetic fields. 

Our analytical description of the atmosphere according to Equations (2)-(4) 
yields of course only a rough approximation to observed spectral lines. In order 
to assess this approximation we compare the synthesized Stokes-/ line profiles 
of the Fe i 6302.5 line for a line of sight in the ambient medium and on the axis 
of the flux concentration, and the profile averaged over the entire box with the 
corresponding spectral line profile from the Jungfraujoch solar atlas (Dclbouille, 
Roland, and Ncvcn, 1973) in Figure 6. The synthetic spectral lines were cal- 
culated without including micro turbulence and hence they are narrower than 
the observed line. The average profile is weaker than the profile emerging from 
the ambient medium in accordance with the observed "line weakening" or "line 
gap" in magnetic elements (Shceley, 1967). Although the Zeeman- broadening 
is strong, we found (by arbitrarily setting the magnetic field to zero for the 
radiation transfer) that the line weakening is mainly due to the strong difference 
in the temperature as a function of optical depth between the magnetic element 
and the ambient medium. This is in agreement with the findings of Shclyag 
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Figure 6. Stokes-/ profiles for the spectral line Fe I 6302.5 A for a line of sight in the ambient 
medium (dashed) , on the axis of the flux concentration (dot-dashed) , and the profile averaged 
over the entire box (dotted) along with an observed profile (solid) from the atlas of Dclbouillo, 
Roland, and Neven (1973). 

et al. (2007) from three-dimensional MHD simulations. The central line of sight 
produces a Stokes /, which is much weaker than the ambient profile and is 
fully split. Indeed, a similar behavior is seen in the observed Stokes-/ profile 
from the fully resolved magnetic element of Lagg et al. (2010). Since our aim is 
to understand the basic principles governing the response of Stokes-X^ profiles 
to the propagation of MHD waves and not to perform a detailed comparison 
with observed Stokes-V^ profiles, we judge the chosen analytical form of the 
atmosphere as good enough and do not strive for a better fit of the synthetic 
with the observed Stokes / by introducing microturbulent and mesoturbulent 
velocities. 

The Stokes-F profiles were computed along the vertical lines of sight by 
integration of the radiative transfer equation for polarized light. If the profiles 
emerging from the top boundary are spatially averaged over the entire width 
of the box, the profiles do not show significant variation with time revealing 
no sign of wave propagation inside the box. This is due to the fact that the 
anti-symmetric flow pattern within the flux sheet will average out to give zero 
net contribution to the Stokes-F variation. This is different in the case when 
the horizontal integration is carried out over only a narrow spatial window 
on either side of the flux sheet axis. In this case, profiles show signatures of 
wave propagation. Therefore, it is necessary to observe at very high spatial 
resolution in order to study the effect of wave propagation in individual flux 
concentrations. The lines of sight that cover only one half of the flux sheet or a 
small part of it give more information about the wave activity in the domain. In 
order to quantitatively study the signatures of wave propagation, we study the 
evolution of the wavelength shift of the central zero crossing of Stokes V (^A^c) 
and the area and amplitude asymmetries, SA and Sa, respectively, according to 
Equations (10)-(12). The zero-crossing shift of Stokes V is defined as, 

6X,c = X.c-Xo, (10) 
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Figure 7. Stokes- profiles of Fc I (a) A 5250.2 A, (b) A 5247.06 A, (c) A 6301.5 A, and (d) 
A 6302.5 A at an elapsed time of 40 s from the vertical lines of sight in a slice ranging from 
X = 410 km to X = 610 km (left of the axis). The magnetic flux sheet has a field strength of 
1000 G on the axis at z = 0. The solid vertical line marks the unshifted central wave-length 
position, the dotted line the zero-crossing position of Stokes V. 



where A^c is the wavelength of the central zero-crossing of the Stokes V profile 
and Ao a reference wavelength, which for the prciscnt purpose is the rest wave- 
length of the spectral line. The asymmetries between the blue and red lobe areas 
of the Stokes-y profiles, and A^, and the amplitudes of the blue and the red 
lobe, ttb and Ur, are defined as, 

M= |t|"|tp (11) 



and 



Sa=\^^^. (12) 
|ab| + l^rl 



5.1. Moderate Field Case 

Figures 7 and 8 show snapshots of Stokes-F profiles of Fe i AA 5250.2, 5247.06, 
6301.5, and 6302.5 A, at time t = 40 s after the start of the simulation of 
the moderate field case (see Figure 4 for reference). Figure 7 shows the Stokes-F 
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Figure 8. Stokes-V profiles of Fe I (a) A 5250.2 A, (b) A 5247.06 A, (c) A 6301.5 A, and (d) 
A 6302.5 A at an elapsed time of 40 s from the vertical lines of sight in a slice ranging from 
X = 670 km to a; = 870 km (right of the axis). The magnetic flux sheet has a field strength of 
1000 G on the axis at z = 0. The solid and dotted vertical lines are as in Figure 7. 



profiles averaged over a horizontal distance from x = 410 km to x = 610 km (left 
of the symmetry axis) and Figure 8 shows the profiles averaged over x — 670 km 
to X ~ 870 km (right of the symmetry axis). The lines Fe i A 5250.2 A and Fe i 
A 5247.06 A belong to the same multiplet of iron, differing only in the effective 
Lande factor, which are 3 and 2, respectively. Hence, the Stokes-^ amplitudes are 
different for the two lines and scale approximately according to the ratio given 
by the Lande factor as 3:2. This can be seen in the plots shown in Figure 7, 
where the amplitude of Fe i A 5247.06 A is lower than that of Fe i A 5250.2 A. 
Similarly, the amplitudes of Fe i A 6301.5 A and Fc i A 6302.5 A scale according 
to the Lande factors of 1.67 and 2.5, respectively. However, the amplitude ratios 
are not strictly according to the Lande factor ratios in the strong-field regime 
because of saturation effects (Stenflo, 1994). 

Apart from the thermodynamic properties of the plasma wherein the spectral 
lines are formed, the height of formation of the spectral lines also depends on the 
line strength, which in turn depends on the excitation potential and the log (gf) 
value. The region of formation spans a wide range and is affected by the presence 
and the height variations of the magnetic field — the latter matters in particular 
for Stokes parameters Q, U, and V. The line depression contribution functions 
for the four Fe i lines are plotted in Figure 9 for the initial model with moderate 
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a) b) 

Figure 9. The line depression contribution functions for the four Fe I hnes of Table 2, 
evaluated for a line of sight located at x = 520 km, in the initial model with 1000 G field 
strength at z = km. (a) Line depression contribution functions for Stokes / in the line core, 
(b) Corresponding contribution functions for Stokes 1/ at A = Ay^. ^ i.e., the wavelength of 
the blue Stokes-V peak. 

field strength and for the hne of sight &l x — 520 km, i.e., in the centre of the 
bundle of lines of sight that was considered in Figure 7. Figure 9a shows the 
Stokes-/ line depression contribution functions in the line cores and Figure 9b 
shows the Stokes-T^ line depression contribution functions at the wavelength 
position of minimal V signals at Av;„i„. For the definition of the Stokes line 
depression contribution function we refer to Grossmann-Doerth, Larsson, and 
Solanki (1988a). Clearly, the maximum contribution to both Stokes / and V 
comes for the doublet Fe i 5250.22 A and Fe i 5247.05 A from higher layers in 
the atmosphere compared to Fe i 6301.50 A and Fe i 6302.50 A. This difference in 
the relative height of formation of these lines is consistent with the theoretical 
models of Khomenko and Collados (2007) and multi-line spectropolarimetric 
observations of Socas-Navarro et al. (2008). 

The asymmetry in the Stokes-^ profiles after 40 s can be clearly seen in the 
Fe I 6301.5 and 6302.5 A lines plotted in Figures 7 and 8. The effect of the 
wave propagation is first sensed by these lines as they are formed lower in the 
atmosphere than the other two lines. The Stokes-F asymmetries as functions 
of time for the above four lines give a clearer picture. Figure 10 shows (a) the 
Stokes-y amplitude asymmetry and (b) the Stokes-F area asymmetry for the 
four Fe i lines as functions of time. The red and blue colours represent the narrow 
bundles of lines of sight on the two sides of the flux-sheet axis. The blue solid 
curves correspond to lines of sight spanning x ~ 410 km to x = 610 km (left) 
and the blue dashed curves for x — 120 km to x — 320 km (far left). Similarly on 
the right side, the red solid curves are for x = 670 km to x = 870 km (right) and 
the red dashed curves are for x = 960 km to a; = 1160 km (far right). In order 
to explain Figure 10, it is essential to look at the velocity field of the simulation, 
which is shown in Figure 4. There, the vertical lines mark the boundaries of the 
line-of-sight bundles. 
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5.1.1. Amplitude Asymmetry 

Initially, there are no velocities inside the domain and hence the Stokes profiles 
are anti-symmetric and therefore 5a is zero. After the start of the simulation, 
velocities start building up in the domain, consequently giving rise to asym- 
metries in the Stokes profiles. Let us now consider the bundle of lines of sight 
on the left side of the axis and the response in 6a for the lines Fe i 6301.5 or 
6302.5 A as shown in Figure 10a. After 20 s, there develops a small downdraft 
within the (blue) bundle of lines of sight in the magnetic region. This leads to 
a red shift of the Stokes-F profile that forms in this region, viz., in the lines 
Fe I 6301.5 and 6302.5 A. The spectral line emerging from the quasi-static, 
partially field-free layer further below is less affected by Doppler shifts. It causes 
an asymmetric illumination of the two flanks of the red-shifted line contribution 
formed further above. This results in the blue lobe of the emerging Stokes- 
V profile be suppressed and consequently the amplitude asymmetry tends to 
become negative. The amplitude asymmetry starts decreasing and reaches a 
minimum at around 50 s, after which time the downdraft moves out of the 
line formation region and a following updraft starts to dominate (see Figure 4), 
making the blue lobe strong again and the amplitude asymmetry to rise. 

On the other hand, the line of sight on the right side of the axis shows an 
inverse time dependence, which comes about because the waves are 180° out of 
phase on the opposite sides of the flux-sheet axis (Section 4.1). There develops 
a strong updraft within the line of sight after w 30 s into the simulation. The 
contribution to the line formation that stems from within this updraft is then 
shifted to the blue relative to the line emerging from the quasi-static layer below. 
This situation suppresses the red lobe of the emerging Stokes-T^ profile, as long as 
the latter is mainly formed within the updraft, making the amplitude asymmetry 
to rise towards positive values. This trend is seen until slightly after 40 s when 
the downdraft of a following wave phase replaces the updraft in these regions. 
The maximum value is reached before the time of minimum value of 6a{t) on 
the left hand side because the wave on the right hand sides of the flux-sheet axis 
is preceding the wave on the; left hand side. However, we notice that the two 
curves for 5a{t) are not symmetric (relative to the time axis) — the (red) curve 
for the right side shows initially even a slight trend towards negative values like 
the (blue) curve for the left side and its amplitude is smaller than that for the 
left side. This asymmetry is due to the fact that the excitation of the flux sheet 
is asymmetric too. The flux sheet continuously moves to the right side until it 
comes to a halt after 12 s and 38.2 km to the right of the initial symmetry axis 
(from Equation (9)). This can be seen in Figure 4, where the blue solid and the 
red solid lines of sight sample more peripheral and more central parts of the flux 
sheet, respectively. We found that the opposite temperature perturbations on 
the opposite sides of the flux-sheet axis has only a minor effect on the behaviour 
of 5a{t). This is also true for the transversal magnetic fleld component. 

In the case of lines formed higher in the atmosphere, like the lines Fe i 5250.2 A 
and 5247.06 A, we see a similar behaviour like described above but with a time 
lag, depending on the arrival time of the perturbations. We see a delay of 9 s 
between the two pairs of lines corresponding to roughly a distance of 58 km. 
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since the sound speed at these heights is around 6.5 kms^^ (see Table 1). This 
distance corresponds to approximately the distance in logr from —1.8 to —2.1 
of the maximum peaks of the Stokes-V^ contribution functions for the two line 
pairs as can be seen from Figure 9b. 

We now consider the lines of sight further away from the flux-sheet axis. 
These are depicted as dashed vertical lines in Figure 4 and the dashed curves 
in Figure 10a show the corresponding asymmetry 6a. The perturbations become 
significant in this region after 50 s only when the first front arrives. On the 
right side of the flux sheet, we see that the velocities are directed upward. This 
suppresses the red lobe making the amplitude asymmetry positive. On the left 
side of the flux sheet, the downward velocities shift the amplitude asymmetry 
towards negative values. We notice, however, that initially there exists the op- 
posite tendency: especially the right side showing negative 6a. This behaviour 
is even more pronounced for the inner bundles of lines of sight on the right 
hand side and corresponding red solid curves of Figure 10a. It indicates that the 
explanation provided above for the origin of the asymmetries is not complete 
even though it can serve as a rough guideline. The origin of this tendency to 
opposite asymmetries is that, prior to the wave affecting the formation of Stokes 
V, it affects the field-free layers below, where Stokes / is already forming. In 
particular, the updraft in the leading wave on the right hand side first causes 
Stokes / to be blue-shifted, which suppresses the blue lobe of Stokes V formed 
higher up in the magnetic region, which is not yet affected by the wave. This 
causes 6a to become negative. Since the velocities are still moderate at this stage, 
the asymmetry remains moderate as well. As the wave moves further up it grows 
in amplitude and enters the magnetic region in the bundles of lines of sight on 
the far right side. This leads to a blue shift of the Stokes-F contribution, while 
the velocities in the field-free layer below are decreasing as the wave moves out of 
this region. Hence, the red lobe of Stokes V gets suppressed. Not until then, 6a 
rises to the expected positive values and these values become substantial because 
the velocities are growing rapidly. 

5.1.2. Area Asymmetry 

The gradients in velocity and magnetic field cause the asymmetry. A rule to 
calculate the sign of the area asymmetry, 5 A, for a purely longitudinal component 
has been provided by Solanki and Pahlkc (1988) (see also Stcincr, 1999): 

A\B{t)\ dz^r<0 M>0, 
dr ■ dr \ > ^ M < 0. 

By convention, v{t) is taken to be positive for flows in the direction of increasing 
optical depth and vice versa, where v is the line-of-sight velocity. In case of a 
flux tube expanding with height, a line of sight along the tube axis will have 
d|i?(T)|/dr > 0. Given this information, together with Equation (13), the correct 
interpretation of the time dependence of 6A would be relatively straightfor- 
ward, except that for a propagating wave the term dv{T)/dT changes sign for 
each half wave. If the line of sight is eccentric, it may traverse the flux-sheet 
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Figure 10. The Stokes-V (a) amplitude asymmetry and (b) area asymmetry for the four Fe I 
Unes hsted in Table 1 as functions of time for the moderate field case with 1000 G. The red 
solid curves represent the bundle of lines of sight on the right side of the axis. The red dashed 
curves represent the bundle of lines of sight on the far right. The blue solid curves are for the 
left slice and blue dashed curves are for the far left bundle. 
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Figure 11. Stokes-V profiles of Fe I (a) A 5250.2 A, (b) A 5247.06 A, (c) A 6301.5 A, and (d) 
A 6302.5 A at an elapsed time of 50 s from the vertical lines of sight in a slice ranging from 
X = 120 km to X = 320 km (far left of the axis). The magnetic flux sheet has a field strength 
of 1600 G on the axis at 2 = 0. The solid and dotted vertical lines are as in Figure 7. 



boundary, where the field drops suddenly with increasing optical depth, so that 
d\B{T)\/dT < 0. In this case, a positive 6A is realized if dw(r)/dr > 0, i.e., if 
there is an accelerating downflow. This is for instance the case when there is no 
velocity inside the flux sheet but a downflow in the field-free surrounding region, 
or in case with an upflow within the flux sheet and no flow in the outside field- 
free region (Grossmann-Doerth, Schiissler, and Solanki, 1988b). Thus, 6A{t) in 
Figure 10b has contributions from gradients in magnetic field and velocity that 
stem from inside the magnetic flux sheet as well as from the boundary of the 
flux sheet. These contributions may have opposite signs. This is also true for the 
contributions from different phases of the wave. In addition, 6 A is an integral 
quantity with contributions over all wavelength of the spectral line. Hence, SA 
forms over a wide height range and therefore includes more than a single wave 
crest or wave trough. This renders the correct interpretation of 5A{t) in this 
case more intricate than it is for Sa{t). This is in particular true for the inner 
bundles of lines of sight. The outer bundles show a behaviour very similar to 
that of Sa{t). 



20 



Stokes ( Lia,a;ii(ist ICS of wave |)roi)agat ion 



8a.1.33e-01 
8A.1.l6e-01 


t . 50s 
Fe 1 5250.22 : 







Sa - 1 .49e-01 
5A.l.18e-01 



1.50s 
Fe I 5247.06 




a) 



b) 



6a = 2.35e-01 
5A . 9.25e-02 


t = 50s 
Fe 16301.50 J 








d) 



Figure 12. Stokcs-V profiles of Fc I (a) A 5250.2 A, (b) A 5247.06 A, (c) A 6301.5 A, and (d) 
A 6302.5 A at an elapsed time of 50 s from the vertical lines of sight in a slice ranging from 
X = 960 km to a; = 1160 km (far right of the axis). The magnetic flux sheet has a field strength 
of 1600 G on the axis at z = 0. The solid and dotted vertical lines are as in Figure 7. 



5.2. Strong Field Case 

The emergent Stokes- y profiles were also computed for a flux sheet with a 
magnetic field strength of 1600 G on the axis at 2: = km. Of the two cases 
shown in Figure 5, we consider in the following only the case with the wide 
excitation region. When spatially averaging these profiles over the entire width 
of the box, they do not show significant variation with time, revealing no sign 
of wave propagation inside the box, similar to the case with a field strength 
of 1000 G. Profiles averaged over smaller slices on cither side of the axis show 
signs of wave propagation, once again emphasizing that lines of sight that are 
placed away from the symmetry axis of the flux sheet yield more information 
about the wave activity. Here, we carry out a similar study of the evolution of 
the Stokes-F asymmetries for the 1600 G case as was done for the moderate 
field case. Figures 11 and 12 show the Stokes-F profiles at time f = 50 s of the 
four Fc I lines under study. Figure 11 shows the Stokes-V profiles averaged over 
a horizontal range from x = 120 km to a; = 320 km (far left bundle of lines of 
sight) and Figure 12 shows the profiles averaged over a range from x = 960 km 
to a; = 1160 km (far right bundle of lines of sight). 
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Figure 13. The Stokes-V (a) amplitude asymmetry and (b) area asymmetry for the four Fe I 
lines as functions of time for the strong field case with 1600 G. The red solid curves represent 
the right side of the flux-shect axis. The red dashed curves represent the bundle on the far 
right. Blue solid curves are for the left bundle and the blue dashed curves are for the far left 
bundle. 
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The Stokes-y amplitude asymmetry and area asymmetry as functions of time 
for these Hnes arc shown in Figure 13. The colour coding is the same as in 
Figure 10. The temperature perturbation and the velocity field for different time 
instances are shown in Figure 5b. There, the colours represent the value of AT 
and the arrows show the velocity vectors at times t = 20, 30, 40, and 50 s. 

Unlike in the case with 1000 G, here we have significant fast, predominantly 
magnetic waves, which get refracted within the flux sheet and convert to fast 
acoustic waves when they encounter the flux-sheet boundary, where they leave 
the flux sheet and enter the field-free domain. This causes the wing-like feature 
in the temperature perturbations of Figure 5b that extends from z = 300 km 
to z = 700 km on both sides of the flux sheet (approximately along the /3 = 1 
contour) at time t = 40 s. The velocities associated with both the fast and the 
slow waves result in the shift and the asymmetries of the Stokes-F proflles. 

In the following, we consider a bundle of lines of sight on the far right side 
of the flux-sheet axis. Here, the front of the fast, predominantly magnetic wave 
starts to become effective at a t w 40 s. The velocities are directed downwards in 
the magnetic region, resulting in a red-shifted Stokes-I^ contribution illuminated 
by light from the unshiftcd absorption formed in the static layer below, which 
suppresses the blue lobe, leading to a negative asymmetry (dashed red curve 
in Figure 13a). At t = 60 s, this downflow has moved out of the line-of-sight 
bundle into the field-free region while within the magnetic region an upflow 
evolves. This upflow, which is due to the following slow mode, gives rise to a 
positive asymmetry leading to a strong positive bump of the red dashed curve 
around t = 70 s. This behaviour is different from the 1000 G case shown in 
Figure 10a wherein the first negative bump due to the fast converting mode is 
less pronounced or missing because of the relatively weak magnetic field. The 
line of sight on the left hand side of the flux-sheet axis shows a similar behaviour 
but of opposite sign. This behaviour can also be seen in the area asymmetry of 
the outer bundles of lines of sight. Thus, 6a and 6 A show a clear signature of 
both the fast and the slow mode. 

5.2.1. Zero-Crossing Shift 

The flows within the magnetic elements can be estimated by the Stokes-F zero- 
crossing shift, given by Equation (10). The 180° out-of-phase flow pattern formed 
on the two sides of the flux-sheet axis creates phase shifted, opposite zero- 
crossing shifts on the two sides. Figure 14 shows the Stokes-F zero-crossing 
shift as a function of time for the moderate and the strong field cases. In case of 
the moderate field and line-of-sight bundles close to the sheet axis, the upfiow 
on the right side and the downflow on the left side of the axis, which stem 
from the respective compressional and rarefactional fronts of the slow (acoustic) 
wave, result in a zero-crossing shift of opposite signs. The curves are shifted in 
time for the four lines and the two line-of-sight bundles due to the difference in 
the formation heights of the lines and due to the time lag of the waves to the 
left and to the right of the flux-sheet axis, respectively. A similar trend in the 
zero-crossing shift can be seen in the outer line-of-sight bundles, where the wave 
fronts arrive later. The flow pattern associated with the fast (magnetic) wave 
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Figure 14. Stokes-V zero-crossing shifts for the four Fe I lines as functions of time for (a) 
the moderate field case with 1000 G and (b) the strong field case with 1600 G. The red solid 
curves represent the slice on the right side of the axis. The red dashed curves represent the 
slice on the far right. Blue solid curves are for the left slice and the blue dashed curves are for 
the far left slice. 
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is a prominent feature in the case with the strong field. It creates a significant 
bump of opposite sign on either side of the flux sheet for the outer line-of-sight 
bundles around t « 40 s, visible in the dashed curves of Figure 14b. 

6. Summary and Conclusions 

This work is an extension of a previous work by Vigeesh, Hasan, and Steiner 
(2009), which focused on the dynamics and the energy transport that occur in 
intense flux tubes as a consequence of an impulsive transversal footpoint motion. 
In the present work, we have constructed flux tubes embedded in the photosphere 
and used the results of our simulation to compute the Stokes profiles that emerge 
from the top of the simulation box in order to study observational signatures of 
wave propagation inside the tubes. 

The nature of the excited modes depends on the value of plasma /3 at the 
place where the driving motion occurs. Depending upon the extent of the region 
of excitation, the energies imparted to the different modes vary. When the 
excitation occurs in a high-/? plasma, we observe that the excited modes are 
a slow acoustic wave and a fast magnetic wave that undergo mode conversion 
and transmission across the f3 = 1 layer. In the case of excitation of a low-/? flux 
sheet over a large enough impact area, most of the energy will go to the fast 
(magnetic) mode. If the area of impact is smaller, then most of the energy goes 
into the slow (acoustic) wave, which is channelled up along the flux tube and 
eventually dissipates by shock formation. Hence, the impact of a large granule 
may impart more energy to the fast (magnetic) mode, but the non-magnetic 
atmosphere gains back this energy in the form of a fast (acoustic) wave due to 
refraction of the fast mode and mode conversion. 

The anti-symmetry in the wave pattern with respect to the flux-sheet axis 
gives rise to distinct observational signatures in Stokes V. While the average 
Stokes-^ profile over the whole domain does not show any significant variation 
with time, clear evidence of the wave phenomena can be detected when looking 
at higher resolved lines of sight on either side of the flux sheet. Stokes- profiles 
become asymmetric, showing opposite temporal behaviour on the two sides of 
the flux-sheet axis. Furthermore, effects of refraction of the fast, predominantly 
magnetic wave in the case of a strong magnetic field are clearly visible in the 
Stokes asymmetry and zero-crossing shift as a function of time. Our results show 
a clear signature of the fast, predominantly magnetic wave in these profiles. 
We come to the conclusion that polarimetric signatures of wave propagation in 
magnetic elements can be observed, provided that the spatial resolution is high 
enough so that magnetic concentrations can be resolved into different regions 
within the flux concentration. Although the simulated Stokes asymmetries would 
be detectable with current polarimetric instruments, observations of considerably 
higher spatial resolution, capable of resolving individual flux concentrations, are 
needed in order to reveal the propagation of waves in individual flux concentra- 
tions and to detect the different modes of the MHD waves. We have highlighted 
the importance of using the Stokes-F asymmetries as a possible diagnostic tool 
to study wave propagation in magnetic elements at disk centre. Observations 
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off disk centre would possibly slightly relax the high requirements of spatial 
resolution because the wave pattern can be expected to be less anti-symmetric 
in this case. However, the interpretation of the polarimetric signals as functions 
of time would become more intricate. The analysis in this work is based on 
photospheric lines under the assumption of local thermodynamic equilibrium 
(LTE). Hence our conclusions are not valid for lines formed in the chromosphere, 
since the LTE approximation is no longer valid in this region. A more realistic 
modeling should be carried out in three spatial dimensions and include NLTE 
effects. 
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